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ABSTRACT 


The  classical  problem  of  celestial  navigation,  in  its 
simplest  form,  is  the  determination  of  an  observer's  longitude 
and  latitude  from  the  altitudes  of  two  identified  stars,  observed 
at  a  known  Greenwich  Mean  Time  on  a  known  date.  A  novel  solution 
of  this  problem  in  closed  analytical  form  is  given  herein.  The 
solution  yields  the  two  possible  positions  of  the  observer  without 
any  prior  knowledge  of  his  position,  without  any  dependence  on 
tables  of  computed  altitudes  and  azimuths  for  an  assumed  position, 
and  without  any  graphical  work.  The  basic  two-fold  ambiguity  is 
resolved  to  yield  a  single,  unique  position  by  repeating  the  calcu¬ 
lation  using  the  altitudes  of  a  third  star  and  of  one  of  the  two 
previous  stars.  The  full  analytical  solution  is  given  as  are 
some  artificial  numerical  examples,  readily  performed  on  a  hand¬ 


held,  programmable  calculator. 
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1.  INTRODUCTION  AND  SUMMARY 


The  classical  problem  of  celestial  navigation,  in  its 
simplest  form,  is  the  determination  of  an  observer's  longitude 
and  latitude  from  the  altitudes  of  two  identified  stars,  observed 
at  a  known  Greenwich  Mean  Time  on  a  known  date.  Numerous  ingenious 
solutions  to  this  problem  by  essentially  trial-and-error  methods 
using  the  concept  of  the  "assumed  position"  have  been  devised  and 
reduced  to  practice  (Bowditch,  1962).  But,  to  the  author's 
knowledge,  no  treatise  on  practical  navigation  draws  attention 
to  the  fact  that  the  basic  problem  can  be  solved  in  closed 
analytical  form  without  reference  to  tables  of  computed  altitudes 
and  azimuths,  graphical  approximations,  etc.  At  least  two 
analytical  solutions  to  this  problem  have  been  published 
(Kotlaric,  1971-72)  (Watkins  and  Janiczek,  1978-79)*  The 
solution  given  herein  appears  to  have  merit  because  of  its 
geometric  simplicity,  its  freedom  from  trigonometric  ambiguities, 
and  its  ready  adaptability  to  a  hand-held,  programmable  calculator. 

Data  from  one  pair  of  star  sights  yield  two  possible  posi¬ 
tions,  usually  widely  separated  and  easily  distinguished.  Data 
from  a  second  pair  of  star  sights  (i.e.,  using  the  altitudes 
a  third  star  and  of  one  of  the  two  previous  stars)  resolve  the 
basic  ambiguity  because  only  one  of  the  two  positions  from  the 
second  calculation  is  identical  to  one  of  the  two  positions  from 
the  first  calculation.  No  prior  knowledge  of  the  approximate 
position  of  the  observer  is  required. 
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2.  STATEMENT  OF  THE  IP OB LEM 

Observed  Data 

hi  ,  h.2  Altitudes  (corrected  for  atmospheric  refrac¬ 
tion,  etc. )  of  stars  1,  2 

t  Greenwich  Mean  Time  of  observation  (on  a 

known  date) 

For  simplicity,  it  is  assumed  that  the  two  observations  are  simul¬ 
taneous. 

Data  Derived  from  Observations  and 
from  the  AIR  ALMANAC 

\\ ,  \2  East  longitudes  of  substellar  points  of 
stars  1,  2 

6i ,  62  Declinations  of  stars  1,  2 
Ci >  £2  Zenith  distances  of  stars  1,  2 

Problem 

To  find  the  Longitudes,  A  (East  +  or  West  -),  and 
Latitudes,  $  (North  +  or  South  -)  of  both  intersections  (P  and  Q) 
of  the  two  small  circles  of  position  on  a  sphere  of  unit  radius 
(Figure  l). 
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3.  AN  ANALYTIC  SOLUTION 

(a)  Each  circle  of  position  is  the  intersection  of  the 
circular  cone  of  half  angle  Q  with  the  unit  sphere  (Figure  2). 

The  center  of  the  circle  has  the  spherical  coordinates  of  the 
star  6  (Figure  3)*  The  equation  of  the  plane  containing  the 
circle  of  position  is  of  the  form 

x  cos  a  +  y  cos  @  +  z  cos  y  -  p  =  0  (l) 

where  a,  3,  and  y  are  the  angles  to  the  X,  Y,  and  Z  axes, 
respectively,  of  the  normal  to  the  plane,  this  normal  being 
through  the  origin  of  coordinates  0,  and  p  being  the  distance 
along  that  normal  from  0  to  the  plane  (Figures  2  and  4).  The 
adopted  XYZ  coordinate  system  is  a  right-handed  one  with  0  at  the 
center  of  the  earth,  with  Z  the  earth’s  polar  axis  (+Z  to  the 
North  Pole),  with  the  XY  plane  the  equatorial  plane,  and  with 
+X  through  the  Greenwich  meridian.  Note  that 

p  «  cos  £ 

and  that 

y  *  90  -  6  . 


Thus 


cos  a  =  cos  \  cos  6 
cos  &  =  sin  \  cos  6 
cos  y  =  sin  6 

After  substituting  the  above  quantities  in  equation  (l),  the 
equation  of  the  plane  is  found  to  be: 

x  (cos\^  cos  6^)  +  y  (sin  cos  6^) 

+  z  sin  6 .  -  cos  C.  =  0, 

i  i 

i  -  1,  2  .  (2: 

(b)  The  two  planes  intersect  in  the  line  PQ. 


Let  a^ 

cos  \ ^  cos 

6. 

i 

b. 

l 

= 

sin  cos 

l 

6i 

c . 

i 

sin  6^ 

pi 

= 

cos  Q± 

The  equation  of  this  line  of  intersection  is  obtained  by 
the  simultaneous  solution  of  equations  (2)  for  i  =  1  and  i  =  2  to 
give: 


-Bz  +C 


-Ey  +  F 
D 


wherein 


A  3  (aiba  -  aabi ) 
B  3  (bsci  -  bica ) 
C  H  (bapi  -  bipa) 
D  3  (aiCa  -  aaci) 
E  3  (bxca  -  baCi ) 
F  3  (caPi  -  cxpa) 


mmmmm 
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Equations  (3)  may  be  regarded  as  parametric  equations  of  the 
line  PQ  with  x  the  parameter  that  specifies  the  position  of  a  point 
on  that  line.  The  other  two  coordinates  of  that  point  are  given  by: 


y 

z 


F  -  Dx 
E 

C  -  Ax 
B 


oo 


(c)  The  two  possible  positions  of  the  observer,  P  and  Q, 
are  the  intersections  of  the  above  line  with  the  unit  sphere; 
that  is,  they  are  the  two  points  for  which 


X2  +  y2  +  z2  =  1  (5) 

(d)  Substitutions  from  Equations  (4)  into  Equation  (5) 
yield  a  quadratic  equation  in  x,  which  has  two  real  roots,  Xp 
and  xA. 

Hr 

(e)  Equations  (4)  are  then  used  to  find  the  other  two 
coordinates  of  P  and  Q: 


p  :  (Xp,  yp,  Zp) 

«  ■  <V  rv  \) 

(f)  Finally,  the  rectangular  coordinates  of  P  and  Q  are 
converted  to  spherical  coordinates  to  yield  the  desired  longitude 
and  latitude  of  each  point,  viz. 

P  5  ^ 


Q  !  V 


Q.E.D 
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4.  ARTIFICIAL  EXAMPLES 


0000  GMT  of  1  September  1975 


(a)  Input  Data  for  A  =  -917532,  $  =  +417662 


Arcturus 


Altair 


\  =  - 125 f 915 
6  =  +  19.317 
C  =  36.704 

X  =  -  427156 
6  =  +  8-799 
c  =  54.382 


Antares 


X  =  -  927581 
6  =  -  26.576 
C  =  68.045 


X££fL  x  =  -  6o.°520 

6  =  +  58.759 
C  =  25.731 


(b)  Calculated  Positions  from  Pairs  of  Star  Sights 


Arcturus  -Altair 

A  =  -  917532 
=  +  4l.66l 

Aq  =  -  957605 
-  2.148 

Arcturus -Antares 

JC  =  -  917532 

=  +  41.662 

Aa  =  -1977841 

=  +  0.156 


Arcturus-Vega 

A  =  -  91?532 
§p  =  +  4l.66l 

A Q  =  -  86?950 
§£  =  +  29.334 

Vega-Antares  =  _  91o552 

$p  =  +  41.662 

A  =  -  42.°  186 

=  +  21.009 

Vega-Altalr  =  _  91?552 

*p  =  +  41.662 

AQ  =  -  55-°550 
=  +  62.295 

Altair-Antares 

A  =  -  91?552 
=  +  41.662 

An  -  -  11?O07 

=  -  37.143 

(c)  Remarks 

It  is  clear  that  one  of  the  two  calculated  positions  in  each 
of  the  six  cases  is  identical  to  the  position  that  was  used  in  the 
construction  of  these  examples.  Ihis  is  the  only  solution  that  is 
common  to  two  or  more  cases. 

The  method  is,  of  course,  applicable  to  observations  of  the 


sun,  moon,  and  planets  as  well  as  to  observations  of  stars 
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CAPTIONS  FOR  FIGURES 

Figure  1.  This  is  a  sketch  of  two  circles  of  position  on  a  unit 

sphere  whose  center  is  the  center  of  the  earth.  The  centers 
1  and  2  of  the  two  circles  of  position  have  the  spherical 
coordinates  \i,  6i  and  \2,  62  of  the  substellar  points  of 
stars  1  and  2,  respectively.  The  line  PQ  is  the  inter¬ 
section  of  the  two  planes  containing  the  respective  circles 
of  position.  The  intersections  of  this  line  with  the  unit 
sphere  yield  the  two  possible  positions  of  the  observer, 

P  and  Q. 

Figure  2.  This  diagram  shows  a  cross-section  of  the  unit  sphere 
that  contains  the  origin  0  and  the  substellar  point  S  of 
an  observed  star.  The  circle  of  position  is  the  inter¬ 
section  with  the  sphere  of  the  circular  cone  whose  axis  is 
OS  and  whose  half  angle  is  £.  The  normal  to  this  plane  OR 
has  length  p. 

Figure  3.  The  east  longitude  X  and  declination  6  (or  latitude) 

of  the  substellar  point  are  shown  relative  to  the  coordinate 
system  in  which  +Z  is  toward  the  north  celestial  pole,  the 
XY  plane  is  the  earth’s  equatorial  plane,  and  the  +X  axis 
lies  in  the  Greenwich  meridian. 

Figure  4.  In  the  same  coordinate  system  as  that  in  Figure  3,  a,  P, 
and  y  are  the  angles  the  line  from  0  to  the  star  makes  with 
the  X,  Y,  and  Z  axes,  respectively. 


